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The cores of ultra-luminous infrared galaxies (ULIRGs) are very dense environments, with a 
high rate of star formation and supernova explosions. They are thought to be sites of cosmic- 
ray acceleration, and are predicted to emit 7 -rays in the GeV to TeV range. So far, no ULIRG 
has been detected in 7 -rays. Arp 220, the closest ULIRG to Earth, has been well studied, and 
detailed models of 7 -ray production in this galaxy are available. They predict a rather hard 7 -ray 
spectrum up to several TeV. Due to its large rate of star formation, high gas density, and its close 
proximity to Earth, Arp 220 is thought to be a very good candidate for observations in very-high- 
energy (VHE; 100 GeV - 100 TeV) 7 -rays. Arp 220 was observed by the VERITAS telescopes 
for more than 30 hours with no significant excess over the cosmic-ray background. The upper 
limits on the VHE 7 -ray flux of Arp 220 derived from these observations are the most sensitive 
limits presented so far and are starting to constrain theoretical models. We also present upper 
limits for the VHE flux from the ULIRG IRAS 17204-0014, the starburst galaxy IC 342, and the 
active galaxy 3C32L 
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1. Starburst Galaxies 

Star formation rates vary from one galaxy to another, as well as for the same galaxy over time. 
Galaxies with a high rate of star formation (compared to other galaxies, or to the same galaxy at 
other times) are called starburst galaxies (SBGs). This period of high star formation generally only 
lasts for a small fraction of the lifetime of the galaxy. Starbursts may be triggered by interactions 
with other galaxies, e.g., mergers. Studying these types of galaxies can provide insight into the 
evolution of galaxies in general as well as our own Milky Way. 

SBGs show strong emission over many wavelengths. For example, there is a well-studied 
correlation between the emission in radio and far infrared (FIR) that is seen in samples of star 
forming galaxies [1], thought to be driven by the star forming rate. A high star formation rate 
means lots of young, massive, stars, which radiate mostly in the UV and have short lifespans 
on the order of tens of millions of years. The UV emission is absorbed by the interstellar dust 
and re-emitted in the FIR. At the end of their lives, many massive stars explode as supernovae. 
Supernova remnants (SNRs) can accelerate electrons and positrons to high speeds, causing them 
to emit synchrotron radiation in radio wavelengths. Thus, galaxies with a high star formation rate 
have high luminosities both in radio and FIR. Galaxies with a FIR luminosity L™ of more than 
10^^ -Lq are also known as ultra-luminous infrared galaxies. 

Similar to the radio-FIR correlation, the existence of a correlation between HE (high energy, 
100 MeV - 100 GeV) and VHE y-ray emission and EIR emission has been suggested as well [1]. 
Again, this correlation is thought to be driven by the star formation rate. Young SNRs can accelerate 
cosmic rays, which produce y-rays by interaction with the ISM. y-rays are either produced by the 
decay of neutral pions, or by cosmic ray electrons (including secondary electrons from pion decays) 
that upscatter IR or radio photons to higher energies. Note that this correlation only strictly holds 
in the calorimetric limit, when the galaxy is dense enough that most cosmic rays interact, rather 
than leave the galaxy [2, 3]. 

Only eight star forming galaxies, four of which are part of the local group, have been detected 
by the Fermi-LKT [2], at energies of 0.1 GeV to 100 GeV. The authors of [2] show that for 
their sample, the HE luminosities were well correlated with the EIR luminosities. None of the 
non-detections of other star forming galaxies seem to contradict this correlation. 

Only two star forming galaxies have been detected in VHE yrays: M82 [4] and NGC 253 [5], 
which is too few to study correlations. However, assuming the VHE-EIR correlation holds, we can 
use the VERITAS measurement of the y-ray flux from M82 to estimate the fluxes from other star 
forming galaxies. Given the luminosity distance Dl, the TeV flux of a star forming galaxy 
can be estimated as 


J FIR r)2 . rTeV 

17 TeV _ ^ ^L,M82 '^M82 

^ ~ n2 ' rFiR 

^M82 

This has been calculated for the three SBGs under study here, see Table |[ Note that these 
estimates are very rough. There have been several efforts to model the cosmic ray acceleration, 
interaction of cosmic rays with the ISM, and production of VHE y-rays in either single galaxies 
[6] or a larger sample of star forming galaxies [1]. Eor example. Table ^ shows the predictions 
for the VHE y-ray flux from Arp 220 by some of those more detailed models. The flux predictions 
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Source 

Dl 

Lg—lOOO^m 

Emin 

Flux (E > Emin) 


[kpc] 

[1O'°-L0] 

[GeV] 

[m-^ s-i] 

M82 [4] 

3.4 

4.6 

700 

3.7 •10-'^ 

IC342 

3.7 

1.4 

500 

~ 1.6-lO-*^ 

Arp 220 

74.7 

140 

500 

~0.4-10-‘^ 

IRAS 17208-0014 

173 

230 

600 

~0.M0-‘^ 


Table 1: Rough predictions for the VHE flux of three starburst galaxies. Luminosity distances Ld and 
FIR luminosities Lg-iooo^im are taken from [2]. The M82 flux as measured by VERITAS [4] was used as 
the reference VHE flux. A spectral index of —2.5 was assumed for all sources, consistent with the M82 
measurements and the expectations for other star forming galaxies. 


range from one order of magnitude below to one order of magnitude above the naive prediction. 
EBL absorption as well as absorption inside the source galaxy are neglected here as they are not 
important for photons with energies below a few TeV for the galaxies studied here. 

1.1 The ULIRG Arp 220 

Arp 220 is the closest ULIRG to Earth, at a distance of about 75 kpc. It is the product of 
a merger of two galaxies and retains two dense nuclei, sites of very high star formation. Due to 
its high density and close proximity to Earth, it is considered a good candidate to search for HE 
and VHE y rays from its star forming regions. Several groups have created detailed models of the 
cosmic-ray acceleration, propagation and interaction in Arp 220 and other galaxies, and used these 
models to predict the VHE /-ray emission [1, 6, 7, 8]. Some of these predictions are collected 
in Table ^ All of the models shown here were tuned to be consistent with the measured radio 
emission. Due to the uncertainties in the astrophysical parameters, there are degeneracies in the 
model parameters that cannot be resolved using this approach. Hence, there is a significant spread 
between the predictions of the different modeling approaches, and also between different par'ameter 
sets for each model. 

The /-ray telescopes H.E.S.S. and MAGIC have looked for VHE /-ray emission from Arp 
220, but neither has detected a significant excess over the background expectations. H.E.S.S. has 
published upper limits on the integral flux (see [9, 10] and Table MAGIC has published upper 
limits on the differential flux [11] which are above even the more optimistic predictions by [6, 8]. 


1.2 The ULIRG IRAS 17208-0014 and the SBGIC342 

IRAS 17208-0014 and IC342 are two star forming galaxies selected as good candidates to 
search for VHE /-ray emission. IRAS 17208-0014 is very bright in EIR and not too far from Earth, 
about twice the distance to Arp 220. IC342 is not nearly as luminous in EIR, but it is very close to 
Earth and was thus also considered as a good candidate for observations. 

The Whipple 10m telescope did not detect any VHE /ray emission from IC 342, leading to an 
upper limit on the integral flux above 600 GeV of 5.5 • 10'^ m^^ s^^ (af 99.9% confidence) [12]. 
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Source 


Emin 

Emax 

Flux (Emin tO Emax) 



[GeV] 

[GeV] 

[m-^ s-i] 

Prediction, 

[1], increased wind, B-field 

500 

— 

* 4 • lO'*^ 

Prediction, 

[ 1 ], standard parameters 

500 

— 

* 1.5- lO'*^ 

Prediction, 

[ 6 ], Blattnig pp cross section 

500 

— 

00 

0 

Prediction, 

[ 6 ], Aharonian pp cross section 

500 

— 

2 • lO'*^ 

Prediction, 

[7], most optimistic 

500 

3000 

2 • 10 - 1 ° 

Prediction, 

[7], most pessimistic 

500 

3000 

5-10-1* 

HESS flux upper limit [9] 

1000 

— 

<3.3-10-‘i 

HESS flux upper limit [10] 

1430 

— 

< 2.5 • 10-1 

VERITAS flux upper limit, this work 

500 

— 

< 2 . 2 - 10 -‘i 


Table 2: Predictions and upper limits on the integral flux of Arp 220 in the VHE range. The two values 
marked with an asterisk (*) were interpolated using the predictions of the fluxes above 300 GeV and above 
1 TeV, assuming a power-law spectrum. Note that the Blattnig parametrization of the proton-proton (pp) 
cross section is now thought to overestimate the VHE 7 -ray flux. 

1.3 The Active Galaxy 3C321 

3C321, the Death Star Galaxy, is a very peculiar system of two galaxies in the process of 
merging. Both have active nuclei, and the larger galaxy has a jet pointing at its companion. 3C321 
has been imaged in many wavelengths from radio to X-ray [13] and there is strong evidence that the 
jet is interacting with the companion galaxy. The observed X-ray emission can be seen as evidence 
of particle acceleration at the interaction sites. Depending on the nature and energy spectrum of 
the accelerated particles and the suiTounding medium, it might be possible for HE or VHE 7 -rays 
to be produced as well. 

There have been no published studies on 7 -ray emission from 3C231 in the HE or VHE range 
so far. 

2. VERITAS observations 

The Very Energetic Radiation Imaging Telescope Array System (VERITAS) [14, 15] is an 
array of four imaging atmospheric Cherenkov telescopes located at the Ered Eawrence Whipple 
Observatory (EEWO) in southern Arizona (31 40N, 110 57W, 1.3km a.s.L). Since its inauguration 
in 2007, it has been upgraded several times. With the current configuration, VERITAS is sensitive 
to air showers induced by 7 -rays with energies between 85 GeV and >30 TeV. Its field of view has 
a diameter of 3.5°, so VERITAS operates in pointing mode: Each night, a number of sources are 
observed for between 15 minutes and a few hours at a time. VERITAS has a single-photon angular 
resolution of less than 0.1° at 1 TeV and an energy resolution of 15 — 25%. Its peak effective area 
is about lO^m^. Due to the large background of cosmic-ray induced showers and the generally low 
fluxes at very high energies, observations over many nights have to be combined to detect weak 
sources. With the current configuration, VERITAS can detect a point-like 7 -ray source with a flux 
of 1% of fhe flux of the Crab Nebula in about 25 hours. 
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Source 

Observing period 

[MJD] 

Livetime 

[h] 

Mean zenith angle 

n 

Arp 220 

55987-56046 

31 

16 

IRAS 17208-0014 

55678 - 56008 

14 

34 

IC342 

55528 - 55572 

3 

37 

3C321 

54503 - 56046 

41 

17 


Table 3: Overview over VERITAS data analyzed for this paper. 


An overview of the data used to obtain the results presented here can be found in Table 
All data were cleaned and analyzed with two independent analysis packages [16], using a Hillas- 
based analysis. The energy reconstruction was done using lookup tables. We applied /-hadron 
separation cuts that were optimized for hard-spectrum sources. All observations were caiTied out 
in wobble mode, with the telescopes pointing 0.5° away from the source, enabling the simultaneous 
estimation of the remaining cosmic-ray background from off-source counts. 

The two objects Arp 220 and 3C321 are separated by about 1° in the sky, and thus can be 
observed simultaneously with VERITAS. However, the sensitivity drops off for larger offsets. The 
analysis of 3C321 uses both data taken pointing at 3C321 (about 10 h in 2008, with VERITAS in 
its original configuration) and pointing at Arp 220 (about 31 h in 2012, after the first upgrade of the 
array). The loss in sensitivity due to pointing away from the source is offset by the larger exposure, 
and the increased sensitivity of the array after the first upgrade. Eor the analysis of Arp 220, only 
data taken pointing at Arp 220 were used, since adding a small data set with poor sensitivity does 
not improve our results. 

3. VERITAS results 

None of the four objects considered for this analysis showed a significant excess above the 
background expectation. Thus, upper limits on the /-ray flux from fhese objecfs have been calcu- 
lafed. Elux limifs and energy fhresholds can be found in Table ^ 

3.1 Upper limit procedure 

We define fhe ON region as a circular region in fhe sky wifh a radius of 0.09°, cenfered on 
fhe source posifion. Even affer /-hadron separation cufs, some level of background confaminafion 
is presenf in fhe ON region. We can esfimafe fhaf level of background by measuring fhe counfs in 
a source-free OEE region wifh fhe same accepfance as fhe ON region. We assume fhaf fhe counfs 
in fhe ON and OEE regions, Aqn and Voff> are disfribufed according fo Poissonian disfribufions 
wifh expecfafion values Xqn =Ns + Nb and Xqff = Here, Ns (Nb) is fhe expected number of 
signal (background) counfs in fhe ON region, and a is fhe ratio befween fhe area of fhe ON and 
OEE regions. 

We can now esfimafe fhe number of signal evenfs presenf in fhe ON region as Ns = Vqn — 
oc-Noff- The /-ray flux would fhen be esfimafed as F = , where Agff is fhe specfral-weighfed 
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Source 

Emin 

Non 

Noff 

a 

Aeff -T 

Flux (E > Emin) 


[GeV] 




[m2 s] 

[m-2 s-1 ] 

Arp 220 

500 

45 

282 

0.17 

8.35 • 10‘^ 

<2.24-10-'^ 

IRAS 17208-0014 

600 

16 

124 

0.17 

4.76 • 10'^ 

< 1.82-lO-*^ 

IC342 

500 

37 

174 

0.18 

0.94 • 10*^ 

<27.6-10-‘^ 

3C321 

500 

52 

363 

0.17 

7.23 • 10*^ 

< 1.90-lO'*^ 


Table 4: 99% confidence level upper limits on the integral flux of the four galaxies studied here, measured 
by VERITAS. A power law spectrum is assumed for all sources, with an index of —3.0 for 3C321 and index 
—2.5 for the three star forming galaxies. 


effective collection area of the array (determined from simulations) and T is the observing time, 
corrected for dead time due to detector readout. 

If there is no significant excess of ON counts over the background expectation, we can instead 
calculate an upper limit on the number of signal events and from that, an upper limit on the 
integral flux via The upper limits shown here are calculated following Rolke et 

al. [17], method 4 (background from sidebands, known efficiency T = T) and reported at 99% 
confidence level. 

3.2 Uncertainties 

There are two sources of uncertainties that have to be considered here: our ignorance of the 
source spectrum and the systematic uncertainty on the VERITAS energy scale/effective area. 

The effective collection area of the VERITAS instrument depends on the source spectrum. Eor 
the three starburst galaxies, it is assumed that the y-ray spectrum should follow the cosmic ray 
spectrum [1, 3, 6, 8]. Eollowing [8], we assume that the spectmm can be modeled as a power law 
with a spectral index of —2.0 to —2.5 at very high energies. Deviations from the power law spec¬ 
trum are expected, e.g., due to absorption in the host galaxy, or by extra-galactic background light, 
but these effects can be neglected below some 10 TeV. We have calculated the upper flux limits 
for several values of the assumed spectral index; here, we quote the weakest limits which were ob¬ 
tained for an assumed spectral index of —2.5. Eor 3C321, there is no published modeVprediction 
of the VHE y-ray spectrum. We calculated upper limits assuming a power-law spectrum with an 
index of —3.0, —3.5, and —4.0. The limits reported here are for an assumed index of —3, which 
gave us the weakest limits of the three test cases. 

There is a 20% systematic uncertainty on the VERITAS energy scale, which translates into a 
30% systematic uncertainty on the integral flux for a spectral index of —2.5, and a 40% uncertainty 
on the flux for a spectral index of —3.0. 

4. Summary and conclusions 

Star forming galaxies, especially SBGs and (U)EIRGs, are expected to emit y-rays in the HE 
and VHE range. Analogous to the radio-EIR correlation, a correlation between EIR emission and 
(V)HE emission has been proposed. Emission in these three regimes is thought to be driven by 
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the star formation rate. The detection of M82 and NGC 235 in VHE by VERITAS and H.E.S.S., 
and of eight star forming galaxies by the Fermi-LAI, supports these expectations. However, in 
order to make sure that our understanding of the underlying processes is correct, it is important to 
investigate a range of galaxies with different properties and star formation rates. 

We have searched for VHE 7 -ray emission from three star forming galaxies with VERITAS: 
the UEIRGs Arp 220 and IRAS 17208-0014, and the SBG IC342. No signal was detected and 
upper limits on the VHE flux from these galaxies were obtained. These upper limits are more 
constraining than previously published observations by other experiments. The non-detection is 
not very surprising, given our sensitivity compared to the predicted fluxes. 

In the case of Arp 220, the VERITAS upper limits are starting to constrain the theoretical 
predictions. We can clearly exclude the more ‘optimistic’ model by [1] at the 99% confidence 
level, even considering the systematic uncertainties of the measurement. However, the upper limits 
presented here are still consistent with the lower range of predictions from this models. There are 
considerable differences between the predictions by different models. Eor example, a very recent 
model [7] puts the predicted emission from Arp 220 about one magnitude lower than previous mod¬ 
els and thus far out of reach for VERITAS. Very deep observations by current generation lACTs 
would be necessary to significantly improve on those flux limits or discover VHE 7 -ray emission 
from Arp 220. Euture experiments with larger effective collection area, such as the planned CTA 
observatory [18], are expected to make definite improvements and either observe the VHE 7 rays, 
or seriously challenge our understanding of the acceleration, propagation and interaction of cosmic 
rays in star forming galaxies. 
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